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Exploring the infrared gluon and ghost propagators using large asymmetric lattices 
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We report on the infrared limit of the quenched lattice Landau gauge gluon propagator computed from large 
asymmetric lattices. In particular, the compatibility of the pure power law infrared solution (^^)^'^ of the Dyson- 
Schwinger equations is investigated and the exponent K is measured. The lattice data favours K ~ 0.52, which 
would imply a vanishing zero momentum gluon propagator as predicted by the Kugo-Ojima confinement mech- 
anism and the Zwanziger horizon condition. Results for the ghost propagator and for the running coupling 
constant are shown. 
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I. INTRODUCTION AND MOTIVATION 

The mechanism of quark and gluon confinement is not fully 
understood yet. The study of the fundamental Green's func- 
tions of Quantum Chromodynamics (QCD), namely the gluon 
and ghost propagators, can help in the understanding of such 
mechanism. Indeed, there are gluon confinement criteria con- 
nected with the behaviour of the propagators at zero momen- 
tum. In particular, the Zwanziger horizon condition implies a 
null zero momentum gluon propagator D{q^), and the Kugo- 
Ojima confinement mechanism requires an infinite zero mo- 
mentum ghost propagator G{q^). The violation of positivity 
for the gluon propagator can also be seen as a signal for con- 
finement 1 29]. 

In QCD, the investigation of its infrared limit has to rely 
on non-perturbative methods like the Dyson-Schwinger equa- 
tions (DSE) and lattice QCD. Both methods have good and 
bad features, namely we can solve analitically the DSE in the 
infrared, but one has to rely on a truncation of an infinite tower 
of equations. On the lattice, one includes all non-perturbative 
physics, but one has to care about finite volume and finite lat- 
tice spacing effects. 

Recent studies of Dyson-Schwinger equations obtained a 
pure power law behaviour for the infrared gluon and ghost 
dressing functions, 
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with K = 0.595 f?l, which implies a vanishing (infinite) gluon 
(ghost) propagator for zero momentum. Studies using func- 
tional renormalization equations Jsl 0] provided bounds in 
the possible values of the infrared exponent 0.52 < K < 
0.595. The infrared exponent obtained from time indepen- 
dent stochastic quantisation Jsl is within these bounds (k — 
0.52145). 

As an infrared (IR) analytical solution of DSE, the pure 
power law is valid only for very low momenta. In figure [2 
the DSE gluon and ghost propagators [6J are compared with 



the corresponding pure power law solution. Note that for the 
gluon propagator the power law is valid only for momenta 
below 200 MeV, and for the ghost the infrared solution is re- 
stricted to even lower momenta. If one wants to use lattice 
QCD to check for a pure power law behaviour, certainly one 
should consider a lattice volume sufficiently large to accomo- 
date a minimum number of points in the region of interest. For 
symmetric lattices this would require a too large volume. A 
cheaper solution is the use of large asymmetric lattices x L,, 
with L, > L,. For example, in flSH [13111 III 

we use a set 

of asymmetric lattices with = 256 and — 8, 10, . . . , 18 
(see 1 10] for the technical details of the simulations). The 
large temporal size of these lattices allow to access to mo- 
menta as low as 48 MeV. Of course, the price to pay are finite 
volume effects caused by the small value of L^. 




FIG. 1 : The DSE solution for the gluon and ghost dressing functions 
ll^ compared with the corresponding pure power laws. 



Recently, some criticism have been raised against the use of 
asymmetric lattices to study the infrared limit of QCD Green's 
functions. In |13], the authors studied the gluon and ghost 
propagators in SU(2) 3-d theory and reported systematic ef- 
fects due to the asymmetry of the lattice, if one compares with 
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a symmetric one. Of course, this comparison cannot be done 
in our case. Also in 1 14] the author reported asymmetry ef- 
fects in the propagators. However, in what concerns the gluon 
propagator, such effects have been already investigated in j^l- 
Furthermore, in |9, 10] it is shown that the approach — > +oo 
is smooth. Indeed, a quick look at the gluon and ghost data in 
fl3ll . again, suggests that the approach +°° is smooth, 
i. e. that the largest symmetric propagators can be obtained 
by extrapolation of the asymmetric ones. In what concerns 
the quantitative results coming from a single large asymmetric 
lattice, all the above studies show that they provide, at least, 
a bound on the infinite volume limit. Here we report on the 
status of our investigations concerning the use of asymmet- 
ric lattices to study the infrared properties of QCD. In what 
concerns the positivity violation of the gluon propagator, our 
results can be seen in L121 . 



II. EXTRACTING THE INFRARED EXPONENT K FROM 
THE GLUON PROPAGATOR 
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FIG. 3: The gluon propagator for all lattices Lj x 256, considering 
only pure temporal momenta. For comparisation, we also show the 
16^ X 48 and 32^ x 64 propagators computed in |2^ . 



In HOT , we have computed the gluon propagator 



(3) 



for SU(3) four-dimensional asymmetric lattices x 256, with 
Lj ~ 8,10,..., 18. To check that the temporal size is large 
enough, 164 16^ x 128 Wilson action gauge configurations 
were generated. In what concerns the gluon propagator for 
time-like momenta, there are no differences between the 16^ x 
128 and 16^ x 256 data (see figure|2}. This gives us confidence 
that the temporal extension of our lattices is sufficiently large. 




q{GeV) 

FIG. 2: The gluon propagator for 16"' x 128 and 16"' x 256 lattices. 



In order to compute the infrared exponent K from the lattice 
data, we considered fits of the smallest temporal momenta of 
the gluon dressing function Zgi^oniq^) to a pure power law 
with and without polinomial corrections. The results are in 
table E] In general, the K values increase with the volume of 
the lattice. So, our K can be read as lower bounds in the infinite 
volume figure K „. 
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0.00 


0.514+ 


12 
15 



0.00 
0.14 
0.18 
0.06 
1.03 
0.00 



TABLE I: K and /d.o.f. from fitting the gluon dressing function 
computed from the different lattices L, x 256. In all fits, the range of 
momenta was chosen such that the fits have one degree of freedom. 
The lowest momentum considered being always the first nonvanish- 
ing momentum. The errors shown are statistical and were computed 
using the bootstrap method, with the number of bootstrap samples 
being about ten times the number of configurations. 



In what concerns the spatial size, it was observed that the 
propagator depends on the spatial size of the lattice (see figure 
|3}- The gluon propagator decreases with the volume for the 
smallest momenta and increases with the volume for higher 
momenta. 



Considering a linear or quadratic dependence on the inverse 
of the volume for the infrared exponent, we can try to extrap- 
olate the figures of table|I] In what concerns the results for the 
pure power law fits, they are not described by these functional 
forms. Using the corrections to the pure power law, one gets 
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values for Koc in the interval [0.51,0.56]; the weighted mean 
value being iCoo = 0.5246(46). 

On the other hand, one can also extrapolate directly the 
gluon propagator, as a function of the inverse of the vol- 
ume, to the infinite volume limit, fitting each timelike momen- 
tum propagator separately. Doing so, we are taking the limit 
Lj OO , assuming a sufficient number of points in the tempo- 
ral direction. Several types of polinomial extrapolations were 
tried, using different sets of lattices, and we conclude that the 
data is better described by quadratic extrapolations of the data 
from the 4th and 5th largest lattices. In figure |3 the two ex- 
trapolations of the gluon propagator are shown. For compari- 
sation, we also include the 16^ x 48 and 32^ x 64 propagators 
computed in i25ll . 

The values of K extracted from these extrapolated propaga- 
tors are K = 0.5215(29), with a j} jd.o.f. = 0.02, using the 
largest 5 lattices in the extrapolation, and K = 0.4979(66), 

jd.o.f . = 0.27 using the largest 4 lattices. Fitting the ex- 
trapolated data to the polinomial corrections to the pure power 
law, one can get higher values for K. Note that the first value 
is on the top of the value obtained from extrapolating directly 
K as a function of the volume. In conclusion, one can claim a 
K e [0.49,0.53], with the lattice data favouring the right hand 
side of the interval. 
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FIG. 4: The bare extrapolated gluon propagator for momenta below 
SOOMeV. The errors on the propagator were computed assuming a 
gaussian error propagation. 



III. ARE DSE FITTING FORMS REPRODUCED BY 
LATTICE DATA? 

In this section we will try to check if the pure temporal lat- 
tice data of the gluon propagator can be described by some 
functional forms that fitted the solution for the gluon propa- 
gator in Landau gauge of the Dyson-Schwinger equations re- 
ported in L15..16.1 . Only the largest lattices will be considered. 



A. The infrared region 

In the IR region, the continuum DSE solution is well de- 
scribed by the following expressions fisll . 
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The first functional form has a branch cut in the IR region. 
The last one has a pole in the IR region. 

In what concerns the lattice propagator, both formulas pro- 
vide good fits up to the maximum of the dressing function — 
see table|n]and figure|5] However, the measured infrared ex- 
ponents are larger than those obtained assuming a pure power 
law for the IR region 1 10.1 . and support an infrared vanishing 
gluon propagator. 
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1.25 
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409: 


H4 
-4 


0.71 


18^ 


x256 
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0.5320+30 


389: 


-6 


1.14 


16^ 


X 128 


570 


0.5100+3? 


416: 


-8 


1.15 


163 


x256 


664 


0.5077+}^ 


409: 


-3 


0.69 


18^ 


x256 


711 


0.5266+2? 


39i: 


-7 


1.09 



TABLE II: IR fits, qmax is the highest momenta considered in the fit. 
qmax and Aqqd are given in MeV. 




q(GeV) 

FIG. 5: The lattice gluon dressing function compared with IR fits. 
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B. The full range of momenta 



The lattice propagator, for the full range of pure temporal 
momenta, was fitted to 



Zfu{q^) = B{q^) {a{q^)) ^ Y= -13/22, 



(6) 



with B[q 



^polexut 



ning coupling were considered lUslfl^ : 



Two different definitions for the run- 
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a(0) + (q^/Alco) X 
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(8) 



The fits using a p{q^) (using po = 1 1) for the running cou- 
pling have a jd.o.f. > 2 for = 16. However, for the 
largest lattice, 18^ x 256, the data is well described; see table 
Uniand figureli 
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X 256 a(0) 


K 




^QCD 


I'/d.o.f. 


B = 
B = 


Zcu, 10.92+J3 


0.5280: 
0.5263: 


h30 
-22 
hl9 
-22 


549+^ 
550+3 


1.61 
1.54 



TABLE III: Fits to all lattice data using ap{q^). 



Note that only when one uses B = Zp„ie the K values agree 
with those from the IR fits. However, AgcD is not compatible 
with the IR values. The fitting form differs from the lattice 
data mainly at the maximum of the dressing function 0.8 
GeV). 

The lattice data adjusts better to (XLN{q^) than to the pre- 
viously considered running coupling. Indeed, looking at the 
fitting results for our largest lattice (see tables Iryl andlVt. the 
values of K and Aqcd are essentially the values obtained in 
the IR study. Moreover, on overall there is good agreement 
between the fitting function and the lattice data, see figure 



ai Tld.o.f. 



16^ X 128 0.5435+^f 364+^ 0.0062+i] 2.44+^ 



16^ X 256 0.5244+^1 374+^ 0.0072+? 2.42+' 



1.82 
1.73 



18^x256 0.5368+^^ 381+^ 0.0067+^ 2.55+^ 1.23 




q(GeV) 



FIG. 6: Compatibility between our lattice data and global fits using 
ap{q^). 
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FIG. 7: Compatibility between our lattice data and global fits using 



A, 



QCD 



ai X^ld.o.f. 



16^ X 128 0.5335+26 3^3+3 O.OO8I+3 2.36+| 
16^ X 256 0.5217+}^ 377+^ 0.0082+;^ 2.36_i 

-17 



s+3 



7+3 

f6 



+ 1 
-1 
+2 



18-' X 256 0.5300+21 388Ij 0.0085+^; 2.46_3 



1.65 
1.61 
1.20 



TABLE V: Fits to all lattice data using UiNiq ) and B = Zpgig. 



TABLE IV: Fits to all lattice data using azA'(? ) and B = 



The running coupling at zero momentum can be measured 
from a.LNi'f') 1 7, 8] and is related with the high momentum 
behaviour of the running coupling, a(0) = (47t/po)fl2 ■ The 
measured a(0) are reported in table I VII Note that the values 



from the 16^ x 128 and 16^ x 256 data agree within one stan- 
dard deviation. Furthermore, a(0) increases with and the 
fitted values reported in table |V1| are around the original DSE 
estimation for a(0) = 2.972 |k2J. 
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163 


X 128 2.79^ 


2 
1 


2.70^ 


-2 
-2 


16^ 


X 256 2.77^ 


1 
1 


2.70^ 


-1 
-1 


18^ 


x256 2.9ll 


2 
3 


2.81^ 


-2 
-3 



TABLE VI: Values of a(0) measured from aiNiq^)- 



IV. GRIBOV COPIES EFFECTS IN THE PROPAGATORS 

The propagators are gauge dependent quantities. To com- 
pute them, one has to choose a gauge. For the Landau gauge 
on the lattice, the procedure consists in maximizing the func- 
tional, 

Fu [g] = Cf £ Re { Tr [g{x)U,{x)g^ {x + fi)] } (9) 

It is well known that this functional has, in general, several 
maxima, the Gribov copies 1 19|. To properly define a non- 
perturbative gauge fixing, one should choose the gauge trans- 
formation g{x) that globally maximizes (|9}. This is a global 
optimization problem and it is not easy to find the global max- 
imum of In most studies of the gluon and ghost propa- 
gators one chooses a local maximum of (|9}, hoping that the 
Gribov copy effects in the propagators are small. 

On the lattice, several studies reported Gribov effects in the 
ghost propagator (see, for example, and the next 

section), but it is generally accepted that Gribov copies do not 
change significantly the gluon propagator. 

However, there are some studies reporting Gribov copy ef- 
fects in the gluon propagator. In lEoll . the problem of the Gri- 
bov copies effects on the gluon propagator was studied, and 
some differences were found. It was claimed a two to three 
O effect due to Gribov copies in the low momentum region. 
Furthermore, it was observed, by sorting the different Gribov 
copies according to the value of the gauge fixing functional, 
that the propagator for the lowest momenta behaves mono- 
tonic ally as a function of < Fu >: for zero momentum, the 
propagator increases with < Fu >, and for some small non- 
zero momenta the propagator decreases with < Fu >. There 
is also a study |21| where the authors claim that the enlarge- 
ment of the gauge orbits allows for Gribov effects in the gluon 
propagator 

In what concerns asymmetric lattices, CEASD gauge fixing 
method i22il was applied to the 16^ x 128 configurations. This 
global optimization method for Landau gauge fixing com- 
bines a genetic algorithm with a local gauge fixing method, 
aiming to find the global maximum of (|9}. The obtained gluon 
propagator was compared with the one computed from a local 
maximum of (|9}; this local maximum is obtained using the 
local gauge fixing method described in lisll . In fig. |S] we can 
see the ratio between the two propagators. 

Although there are many momenta for which the ratio is 
not compatible with one, we cannot conclude in favour of any 
systematic effect. 




ht 



FIG. 8: Gribov copy effects in the gluon propagator computed from 
16^ X 128 configurations. 



V. GHOST PROPAGATOR AND RUNNING COUPLING 
CONSTANT 

A. Ghost propagator 

For our smallest lattices (10^ x 256, 12^ x 256, 16^ x 128), 
we have computed the ghost propagator LlLl . 

G"^{q) = -d"''G{q^) (10) 

using both a point source method i23ll (we averaged over 7 
different point sources to get a better statistics), and a plane- 
wave source 1 17|. In both cases we used the preconditioned 
conjugate algorithm, as described in fl^l . The advantage of 
using a plane-wave source being that the statistical accuracy 
is much better, but we can only obtain one momentum com- 
ponent at a time. Using a point source method one can get all 
the momenta in once, but with larger statistical errors. 

In figures |9] and ^| it is shown the ghost dressing func- 
tion Zghost{q^) for 12^ X 256 and 16-' x 128 computed with 
both methods. For the 12^ x 256, the point source data is not 
smooth for a large range of momenta. This is true also for 
10^ X 256 lattice , but this effect is significantly reduced for 
16-3 x 128. In what concerns finite volume effects, we see 
differences, in the infrared, between pure temporal and pure 
spatial data, as in the gluon case. 

In figure one can see the ghost dressing function only 
for the plane-wave data. As in the gluon case, we are able to 
evaluate the effect of Gribov copies on the lattice 16^ x 128 
by considering different gauge fixing methods. One can see 
clear effects of Gribov copies, as expected from other stud- 
ies fni [T^l . Also, we see finite volume effects if one com- 
pares propagators computed from lattices with different spa- 
tial sizes. 

In what concerns the infrared region, we were unable to 
fit a pure power law, even considering polinomial corrections 
i\ ill . This negative result can be due either to the finite volume 
effects caused by the small spatial extension of our lattices, or 
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FIG. 9: Bare ghost dressing function for a 12 x 256 lattice; "p2p" 
("pws") stands for the ghost components computed using a point 
(plane wave) source. 



FIG. 11: Bare ghost dressing function for all lattices. The data were 
computed using a plane wave method. 
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FIG. 10: Bare ghost dressing function for a 16 x 128 lattice. 



to an insufficient number of points in the infrared region — 
note that for the ghost propagator, the pure power law lacks 
validity well below 200 MeV (see figure^. 



B. Running coupling constant 



and results from lattice simulations fin H m m, shows 

a decreasing coupling constant for small momenta. Using an 
asymmetric lattice allow us to study smaller momenta having 
in mind to provide, at least, a hint to this puzzle. 

Again, our lattice data shows finite volume effects, if one 
compares pure temporal and pure spatial momenta, see fig- 
ure^] Comparing the results for all available lattices (plane- 
wave source), see figure^J we can see, as in the ghost case, 
finite volume effects, and clear Gribov copies effects. 
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Using the gluon and ghost dressing functions, one can de- 
fine a running coupling constant as 



(11) 



The DSE infrared analysis predicts a running coupling con- 
stant at zero momentum different from zero, as{Q) = 2.972 
01 • On the other hand, the DSE solution on a torus i24ll . 



FIG. 12: Running coupling constant for a 16 x 128 lattice. 

In what concerns the infrared behaviour, we tried to fit the 
lowest momenta to a pure power law, {q^)^°-- We concluded 
that this power law is only compatible with the data from 
16^ X 128 lattice, gauge fixed with CEASD method, giving 
Ka ^ 0.688, with X^/d.o.f. - 0.011. The reader should be 
aware that it is also possible, in some cases, to fit the infrared 
data to a(0)(l +aq^ + . . .) and get a a(0) ^ 0. Therefore, we 
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VI. CONCLUSIONS 
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FIG. 13: Running coupling constant for all lattices. The data were 
computed using a plane wave method. 



can not give a definitive answer about the behaviour of the run- 
ning coupling constant for q = Q. Note, however, that 0,s{q^) 
for the smallest momenta, seems to increase as a function of 
the volume. 



In this work asymmetric lattices were used to study the in- 
frared behaviour of QCD propagators. 

In what concerns the gluon propagator, the lattice data is 
well described by a pure power law {q^)^^ for momenta below 
150 MeV. The value of the infrared exponent increases with 
the volume, so our k's can be read as lower bounds of the 
infinite volume figure. Extrapolating the gluon propagator to 
the infinite volume limit, assuming a sufficient temporal size 
of our lattices, one gets K G [0.49, 0.53] . Unfortunately, we can 
not give a definitive answer about the behaviour of the gluon 
propagator at zero momentum. Note, however, that the lattice 
data favours the right hand side of the given interval, and that 
using other fitting forms and a larger range of momenta, one 
always get K > 0.5. 

In what concerns the ghost propagator, we were not able to 
extract a pure power law on the available results. We observed 
finite volume effects and clear Gribov copies effects on this 
propagator, in agreement with other studies. 

Finally, we observed a decreasing running coupling for low 
momenta, in agreement with previous simulations. 
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